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Abstract 



We have performed relativistic many-body calculations of the hyperfine inter- 
action in the 6s and 7s states of Cs, including the off-diagonal matrix element. 
The calculations were used to determine the accuracy of the semi-empirical 
formula for the electromagnetic transition amplitude (6s | Ml \ 7s) induced by 
the hyperfine interaction. We have found that even though the contribution 
of the many-body effects into the matrix elements is ve ry large, the square 
root formula (6s\Hhf s \7s) = J (6s|i// l j s |6s)(7s|i^/ l j s |7s) remains valid to the 
accuracy of a fraction of 10~ 3 . The result for the Ml-amplitude is used in the 
interpretation of the parity-violation measurement in the 6s — 7s transition 
in Cs which claims a possible deviation from the Standard model. 



PACS: 32.80.Ys,31.15.Ar,32.10.Fn 



1 



Typeset using REVT^X 



I. INTRODUCTION 



Recent progress in highly accurate measurements of parity non- conservation (PNC) in 
atoms has got to the point where new physics beyond the Standard Model of elementary 
particles can be studied. The latest analysis of the most precise measurements of the 
PNC in cesium [fj suggests that the value of the weak charge of the 133 Cs nucleus may differ 
from the prediction of the Standard Model. In that experiment ||2j the ratio of the PNC El 
amplitude to the tensor polarizability (3 for the 75*1/2 — 661/2 transition was measured with 
0.35% accuracy. The measured value can be written in the form 

kpNC Qw /-, \ 

(3 N ' 1 j 

where kp^c is the electron matrix element of the electric dipole transition induced by the 
weak interaction between 75i/ 2 and 6S1/2 states of 133 Cs, Qw is the weak nuclear charge 
and N is the number of neutrons. To interpret the measurements in terms of the weak 
nuclear charge one needs to know kp^c and (3. The value of kp^c can be obtained from 
atomic calculations only. Bennett and Wieman [|IJ] used the value /cpnc = 0.9065(36)ieao 
which is the average of our result /cpnc = 0.908(9)ieao obtained in 1989 and the result 
of the Notre-Dame group &pnc — 0.905(9)iea [|J obtained in 1990. Note that Bennett and 
Wieman assumed 0.4% accuracy of the calculations contrary to the 1% accuracy claimed in 
both calculations. This assumption was based on the comparison of the calculated atomic 
quantities relevant to the PNC amplitude (electromagnetic transition amplitudes between 
lower s and p states and hyperfine structure intervals of these states) with the latest very 
accurate measurements which resolved major discrepancies between theory and experiment 
in favor of theory. 

The most precise value of (3,(3 = 27.024(43) (67)ag, was obtained in Ref. Q from the 
measurements of the ratio Ml^/ s //3 where Mlhfs is the Ml transition amplitude between 
the states 6S and 7S induced by the hyperfine structure (hfs) interaction. Semiempirical 
formula for the Ml h f s amplitude derived in Refs. 0-0] was used in the analysis: 



Ml hfs 



c 



^5^1(^-^)1.0024 (2) 



Here Aq s and Aj s are the hfs constants of the 6s and 7s states of Cs, gs = 2.0025, gi = 
—0.0004, the coefficient 1.0024 was introduced to account for the many-body effects. This 
gives Ml hfs = |^|0.8094(20) x 10" 5 §g. 

Values (3 = 27.024(43) (67)a 3 and k PNC = 0.9065(36)iea and measurements of (|TJ) § 
lead to the value of the weak charge of 133 Cs Qw = —72.06(28) (34) which differs from the 
prediction of the Standard Model Q w = -73.20(13) [| by 2.5cr. 

From the point of view of accurate atomic calculations, there are two major questions in 
the analysis above which should be considered. The first is whether the actual accuracy of 
the PNC calculations is really 0.4%. The second is whether the semi-empirical formula (0) 
is accurate. In the present paper we address the second question, leaving the first one for 
later work. 
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II. PRELIMINARY ANALYSIS 



Mlhf s amplitude appears due to mixing of the 6s and 7s states by the hfs interaction, 
Ml hfs = {6s ' F J Hh i 7s ' F) (7s,F\Ml\7s,F>) 



7/ . (6s, F'\H h f s \7s, F') 



+ (6s,F|Ml|6s,F') N ' 1 rajA ' ' ' . (3) 

E-js — tiQ S 

Two major assumptions have been made to arrive at (0) from (|3]). First, the non-relativistic 
expression for the operator of the Ml transition was used: 

Ml = -\fi B \g (L + 2S). (4) 

Second, the square root formula is assumed to be valid 

(6s\H hfs \7s) = yJ(6s\H hf8 \6s)(7s\H hf8 \7s). (5) 

The accuracy of both of these assumptions needs to be examined. The situation is clear with 
the relativistic corrections to the Ml operator (|4]). According to the estimations of Bouchiat 
and Piketty || the relativistic effects modify the amplitudes (6s | Ml | 6s) and (7s | Ml \ 7s) at 
only the 10~ 4 level. This is in line with the many-body calculations of the relativistic effects 
in (/-factors and Ml-transition amplitudes for Cs and other alkaline atoms in our early works 
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The situation with the square root formula (|5|) is less clear. In their pioneering work 
Bouchiat and Piketty estimated the first order core polarization corrections to it and 
introduced the correction factor 1.0017. In a later paper by Bouchiat and Guena pfl this 
factor was assumed to be 1.0024 (see also formula (|2|)). The accuracy of the estimation of 
the many body correction was assumed to be approximately equal to the correction itself 
(~ 0.002) HdJ. In these works there were no accurate calculations of other many body 
contributions to the hfs beyond the first order core polarization corrections. However, it is 
known that these contributions can be up to 20% of the hyperfine structure (see below). 
The applicability of eq. (|J) in this situation is not obvious. 

The accurate relativistic many-body calculations of the off-diagonal hfs matrix element 
(H) were recently performed by the Notre Dame group [|ll[j. The accuracy of the calculations 



was about 1% and agreement with formula (|5|) within this accuracy was achieved. Note that 
the theoretical accuracy for the diagonal hfs matrix elements is also about 1% (see Refs. 
|i"2l , |l3H and this article). This accuracy is not sufficient to find an accurate value of (3 to add 
anything new to the result of the cesium PNC experiment published in 0. 

However, we believe that the validity of the square root formula ([5]) can be demonstrated 
to much higher accuracy than the absolute theoretical accuracy of the hfs calculations (here 
we agree with [0-0]). We suggest that the following combination of matrix elements be 
calculated 

R = {Qs\H hfs \7s} i 

f(6s\H hfs \6s)(7s\H hfs \7s) 
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where all hfs matrix elements are calculated in the same approximation. The value of R 
can be calculated with very high accuracy because uncertainties in different matrix elements 
cancel each other almost exactly We will demonstrate that inclusion of different many body 
and relativistic effects leave the formula 



(6s\H hfs \7s) = ^{6s\H hfs \6s)(7s\H hfs \7s) (7) 

valid to very high accuracy, so that the value of R (§) remains very small. 

Let us start from the analytical estimates of different contributions to R in First 
note that in the single-electron approximation formula (^) is exact if the wave functions of 
the 6s and 7s states are proportional 

= Bi/j 7s (8) 

(a is Bohr radius) on short distances from the nucleus, r < a /Z. Dirac equations for the 
states 6s and 7s differ by the energy only. Therefore, their solutions on short distances 
where the difference in energies is small compared to the potential, differ by normalization 
only. One can say that is valid if 

AE/\V\<£l, (9) 

where AE = 0.08445 a.u. is the energy difference between the 6s and 7s states of Cs, 
V is the atomic potential. The Hamiltonian of the hfs interaction Hhf s is proportional to 
1/r 3 and the main contribution to its matrix elements comes from the distances r < clq/Z. 
Substitution of V = Ze 2 /r, r = a /Z and Z = 55 into ([J) gives 

AE , 

— ^3xl0~ 5 . (fO) 

Note that for s-waves the correction can be even smaller. Indeed, in the non-relativistic 
approximation s-wave hfs is proportional to 5(r). Thus, the typical distances r ~ h/ (m e c) = 
aa , where a = 1/137. 

Let us now consider the many-body effects. It is convenient to do this using the many 
body perturbation theory in the residual Coulomb interaction U, U = H — Hhf- Here H is 
the exact Hamiltonian of the atom and Hhf is the Hartree-Fock Hamiltonian. We generate 
the complete zero-approximation set of the eigenvalues, wave functions and Green's functions 
using the Hartree-Fock Hamiltonian. The small parameter of this many-body perturbation 
theory is the ratio of the non-diagonal matrix element of the residual interaction U to 
the large energy denominator for excitation of the electron from the closed electron shell 
(electron core), e.g. 5p -electron: U/E 5p ~ 10~ 2 . 

The perturbative (correlation) corrections to the hfs matrix element can be divided into 
two classes: the self-energy corrections and the vertex corrections. The former can be 
included into eq. @ through the redefinition of the single electron wave functions while the 
latter are included through the redefinition of the Hhf s operator. 

Self-energy corrections dominate in the hfs of alkaline atoms (see, e.g. [0). The major 
contribution is due to the correlations between an external electron and core electrons. We 
include them by using so called Brueckner orbitals instead of the Hartree-Fock orbitals as the 
single-electron wave functions in eq. (0). The Brueckner orbitals are obtained by introducing 
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an additional operator £ into the Hartree-Fock equations for the external electron and 
solving the Dyson-type equation (Hhf + ^{E) — E)ip = 0. The £ is an energy- dependent 



non-local operator which is also called the "correlation potential" [ll^llql . For the calculation 
of £ see the next section. The Brueckner type correlation correction constitutes 20% of the 
hfs of 6s and 7s states of Cs. However, if we neglect the dependence of £ on energy, 
the estimation ( P~0| ) is still valid. It follows from the calculations that d'E/dE ~ 1% for 
E ~ Eq s ,E 7s (it is suppressed by the parameter AE/E 5p ). This leaves condition (|10"D 
practically unchanged. 

Dominating vertex corrections to the hfs matrix element are due to the effect of core po- 
larization by the nuclear dipole magnetic field. Since the core states change in the magnetic 
field, the Hartree-Fock potential V created by the core electrons as well as the correlation 
potential £ also change. The effect of this change on the hfs can be accounted for by 



redefining the operator of the hfs interaction [12[] : 

H' hfs = H hfs + 5V + 5E. (11) 
The correction to the hfs caused by 5V is often called the RPA-type |16|] correction, while an- 



other correction associated with 5£ is the non-Brueckner correlation correction or structural 
radiation | 15fl . These corrections are more likely to cause deviation from the square root 



formula since they are localized on larger distances up to the core radius. Note, however, 
that in the case of the hfs interaction 5V is completely due to the Hartree-Fock exchange 
potential. There is no change to the Hartree-Fock direct potential since magnetic field does 
not change electron density in the first order of perturbation theory. This means that 5V 
vanishes exponentially outside the core. Inside the core, at r ~ a , AE/V ~ 0.01 and the 
6s and 7s orbitals are still proportional. Note that the potential V at these distances may 
be estimated as V ~ —Z e jfe 2 /r, where Z e ff ~ 5. Since the contribution of SV is about 10% 
we come to the estimate 10~ 3 for the error of the square root formula. 

There is one more reason why the square root formula is accurate. The expression for 
R (^) is symmetric with respect to the energies E 6s and E 7s . Therefore, its decomposition 
over AE (AE = E es - E 7s ) starts from AE 2 : 

R = a(AE) 2 + b(AE) A + (12) 

Since all linear in AE terms are canceled out one can say that the error should be smaller 



than in the estimates above. The dimensionless parameter for (12) is 

(AE/E 5p ) 2 ~ lO" 2 , 

where E 5p « 0.84 a.u. is the core excitation energy. Since the term a(AE) 2 arises due to SV 
and ST, which contribute about 10% and 1%, respectively, into the hfs, the total deviation 
from the square root formula caused by the RPA and non-Brueckner corrections should be 
smaller than 10" 1 x 10" 2 = 10" 3 . We may add that the contribution of 5T, to the hfs is much 
smaller than the contribution of £ since <5£ has an additional suppression by the parameter 
(AE/E 5p ) |§. 

There are also contributions to the self-energy and vertex due to the radiative corrections. 
We have not considered these contributions in our calculations. However, they come from 
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the very short distances r < h/m e c = aa^ and should not cause any significant deviation 
from the square root formula. 

Finally, let us estimate contributions to Mly s which cannot be presented in the form of 
eq. (H). Let us use the basis of the exact atomic eigenstates and treat H h f s as a perturbation 
in this basis. The result can be presented in the form 

s-^ (6s, FlHhtJa, F) , . ,~ 7> 
Ml hfs = J2 ~ F - f - (az,F\Ml\7s, F ) 



-6..s 



(p,F'\H hft \7s,F') 



+ J2(Qs,F\Ml\{3,F y H > ' ' (13) 

Here |6s), |7s), |a), |/3) are the eigenstates which include all possible configuration mixing, 
| a) and \(3) may contain an arbitrary number of pairs of excited electrons and holes in 
the electron core. All non-diagonal matrix elements of the Ml operator vanish in the 
non-relativistic limit. Moreover, it was demonstrated in our work |J that the dominant 
contribution appears only in the second order in the spin-orbit interaction and in the 
first order in configuration mixing, i.e. non-diagonal Ml matrix elements are of the or- 
der Ml ~ (Za) 4 Qi n /E ~ 10~ 4 — 10 _5 |/xb|, where Qi n is the non-diagonal matrix ele- 
ment of the Coulomb interaction corresponding to an excitation of a core electron and 
E is the energy of this excitation. Indeed, the operator of the magnetic moment is 
Ml = /Ub(L + 2S) = /ie(2J — L) (relativistic correction to this expression ~ 10~ 5 ). Elec- 
tron wave functions are the eigenfunctions of the total electron angular momentum J. There- 
fore, J does not give any non-diagonal matrix elements. On the other hand, the matrix ele- 
ment (a, J = l/2|L|/3, J = 1/2) requires spin-orbit interaction both in the bra (a, J = 1/2 1 
and ket \/3,J= 1/2) vectors, since in the non-relativistic limit they correspond to the total 
orbital angular momentum L = 0, i.e. L|a) = L|/3) = 0. Thus, we need the second order in 
spin-orbit interaction. Note that the non-diagonal in angular momentum L matrix elements 
of the hyperfine interaction like (Si/2\Hhfs\d3/2) do not help since in this case both the hfs 
matrix element and Ml matrix element (§1/2! Ml 1^3/2) are very small. 



Non-diagonal matrix elements of Ml were calculated in Refs. the value 

(6s|Ml|7s) ~ 0.4 x 10~ 4 |/is| was measured in Refs. |17|,[5|,[l8|. Thus, each term with the 
non-diagonal Ml matrix element in eq. ( |T3"D is suppressed by a factor of 10~ 4 — 10~ 5 . 
Therefore, we may safely assume that the correction to the diagonal Ml contribution ([|) 
does not exceed 1CT 3 . 

We should note that it may not be easy to come to this conclusion using perturbation 
theory in the Dirac basis of electron orbitals (jj scheme) (Dirac basis was used in Ref. 0). 
In this basis the small result must appear due to strong cancellations between different terms 
in the sum over intermediate states. 



III. MANY-BODY CALCULATIONS 

To test the validity of the square root formula (^) we performed accurate many-body 
relativistic calculations of the off-diagonal and diagonal hfs matrix elements. Detailed dis- 



cussion of the accurate hfs calculations can be found elsewhere |T2|]. Here we repeat the 
main points emphasizing the role of different many-body effects. 
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We start calculations from the relativistic Hartree-Fock (RHF) method in the V N ~ l ap- 
proximation (calculations for the external electron are carried out in the frozen self-consistent 
field of the core ). The core polarization is calculated using the Hartree-Fock equations in 
an external field |14[]. It is equivalent to the well-known random-phase approximation with 
exchange method (see, e.g. |n|). The many-body effects such as the Brueckner-type cor- 



relations, and the structural radiation are included by means of the correlation potential 
method [IS|. As it was pointed out in the previous section, the Brueckner-type correlation 



corrections are included by solving the Dyson-type equation for the states of the external 
electron 

(H HF + £-E)4> = 0. (14) 

Correlation potential E accounts for the correlation between an external electron and core 
electrons. We use many body perturbation theory and the Feynman diagram technique 



to calculate £ [p!2[|19[ . The perturbation expansion of £ in the residual Coulomb interac- 
tion starts from the second order. The corresponding diagrams are presented on Fig. |T|. 
We include both the second order diagrams and three dominating classes of higher order 
correlations: 

1. Screening of the Coulomb interaction between an external electron and core electrons 
by other core electrons. This is a collective phenomenon and the corresponding chain 
of diagrams is enhanced by a factor approximately equal to the number of electrons 
in the external closed subshell (the 5p electrons on Cs). We stress that our approach 
takes into account screening diagrams with double, triple and higher core electron 
excitations in contrast to the popular coupled cluster method where only double and 
selected triple excitations are considered (see, e.g. The effect of screening is 
taken into account in all orders by summation of the corresponding chain of diagrams 
which in the Feynman digram technique form a matrix geometrical progression. 

2. Hole-particle interaction in the core polarization operator. This effect is enhanced 
by the large zero-multipolarity diagonal matrix elements of the Coulomb interaction. 
We take it into account by amending the direct Hartree-Fock potential in which the 
polarization operator is calculated. 

3. Iterations of the self-energy operator (£). This chain of diagrams describes the non- 
linear effects of the correlation potential and is enhanced by the small denominator, 
which is the excitation energy of an external electron (in comparison with the excita- 
tion energy of a core electron). The iterations of £ are included by solving equation 

(ED- 



Substituting the Brueckner orbitals into eq. (|6|) accounts for the dominating correlation 
corrections to the hfs. Corresponding diagrams are presented on Fig. These corrections 
constitutes 23% of the hfs of the 6s state of Cs and 12% of the hfs of the 7s state of Cs. 

To take into account the core polarization effect we self-consistently solve the Hartree- 
Fock equation for the core states in the nuclear magnetic field. The details are presented in 



Ref. flU}] . When all corrections 5ip n to core states caused by the magnetic field are found, 



they are used to calculate the correction 5V to the Hartree-Fock potential. Then core 
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polarization is included into the single-electron matrix element (a\Hhf s \b) between valence 
states | a) and \b) by redefining the operator of the hfs interaction H' h f B = H h f s + SV. This 
corresponds to the summation of the infinite series of the RPA-type of diagrams presented 
on Fig. H The RPA-type core polarization contribution to the hfs of the 6s and 7s states 
of Cs is about 15%. 

Core polarization also leads to the change of E. Corresponding contributions to the hfs 
matrix element (a|5E|6) are often called structural radiation. Second order diagrams for 
structure radiation are presented on Fig. f|. We use direct summation over the complete set 
of single-electron states to calculate these diagrams. 

There is also a contribution to the hfs due to the change of normalization of the wave 
function caused by E. This contribution can be written in a form 

A norm = l -(a\H hfs \b){{a\dt/dE\a) + (b\dt/dE\b)). (15) 

The combined contribution of the structural radiation and renormalization into the hfs of 
the 6s and 7s states of Cs are 1.5% and 0.6% respectively. 

When all dominating higher-order correlations are included into the calculation of the 
Brueckner orbitals for the 6s and 7s states of cesium, the accuracy for the calculated energies 
of these states is very high and constitutes about 0.1%. However, we introduced fitting 
parameters to re-scale E to fit the energies exactly. This procedure allows us to effectively 
include some omitted higher order correlations and to test the sensitivity of the hfs matrix 
elements on the value of E. 

The results for the hfs are presented in Table 0J In this table h = Hhf s , the matrix 
elements of 5V are RPA-type corrections, the matrix elements of 5E are structural radiation 
(including renormalization flT5])), matrix elements with ipBr include Brueckner-type corre- 



lation corrections. One can see that the correction to the square root formula due to the 
considered many-body effects does not exceed 4.4 x 10~ 4 . When all dominating many-body 
effects are taken into account the accuracy of the calculated hfs constants compared to ex- 
periment is about 1%. However the square root formula is still valid to the accuracy of about 
10~ 4 . The most likely cause of the remaining discrepancy with experiment is higher-order 
correlation corrections not included in our calculations. These corrections are localized on 
the radius of the core and due to the fact that these corrections are very small (~ 1% of 
the experimental hfs) it is extremely unlikely that they can break the square root formula. 
The same may be said about the very small radiative and Breit corrections. Note that our 
final results for the diagonal hfs matrix elements for the 6s and 7s states are in very good 
agreement with the calculations of the Notre-Dame group |13| . 



It follows from the above that the correction to the square root formula is about an 
order of magnitude smaller than the estimations of Bouchiat and Guena |7|. However, 
there is no formal disagreement between the results, since Bouchiat and Guena estimated 
the uncertainty of their result to be equal to the correction itself. We believe that for the 
analysis of the PNC experiment it is safer to assume no correction to the square root formula. 
This slightly changes the numbers. The Mlhf s amplitude, tensor polarizability (3 and weak 
charge of the 133 Cs nucleus become 

Ml hfs = |^|0.8074(8) x 10" 5 , 
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(3 = 26.957(43) (27)ag, 
Qw = -71.88(28)(29). 



(16) 



To stress the importance of the result here we used an estimate of the theoretical accuracy 
0.4% in the value of kpNc- Our result for Ml h f s is in very good agreement with the 
result of Derevianko et al |TT| 

Ml hfs = |^|0.8070(73) x 10~ 5 , (17) 

but has the better accuracy. The weak nuclear charge Qw in (0) represents even larger 
deviation from the Standard Model value Qw = —73.20(13) || than the result presented 
by Bennett and Wieman |IJ. The deviation is 2.9<r if 0.4% accuracy of calculations of the 
kpNC is assumed. Note that even if 1% accuracy is assumed for the calculated value of 
kpNc as h was claimed in both theoretical works then there is still 1.5a deviation from 
the Standard Model. However, we would like to stress once more that before making any 
conclusions about agreement or disagreement with the Standard Model the question about 
the accuracy of the atomic calculations of the PNC electronic matrix element kp^c (see (|l|)) 
should be carefully re-analyzed. 



9 



REFERENCES 



[1] S. C. Bennett and C. E. Wieman, Phys. Rev. Lett. 82, 2484 (1999); 82, 4153(E) (1999); 

83, 889(E) (1999). 
[2] C. S. Wood et al, Science 275, 1759 (1997). 

[3] V. A. Dzuba, V. V. Flambaum, and O. P. Sushkov, Phys. Lett. A 141, 147 (1989). 
[4] S. A. Blundell, W. R. Johnson, and J. Sapirstein, Phys. Rev. Lett. 65, 1411 (1990). 
[5] M.-A. Bouchiat, J. Guena, and L. Pottier, J. Phys. Lett. (Paris) 45, L-61 (1984). 
[6] C. Bouchiat and C. A. Piketty, J. Phys. (Paris) 49, 1851 (1988). 
[7] M.-A. Bouchiat and J. Guena, J. Phys. (Paris) 49, 2037 (1988). 

[8] W. J. Marciano and J. L. Rosner, Phys. Rev. Lett. 65, 2963 (1990); 68, 898(E) (1992). 
[9] V. V. Flambaum, I. B. Khriplovich, and O. P. Sushkov, Phys. Lett. 67A, 177 (1978); 

Zh. Eksp. Teor. Fiz. (Sov. Phys. JETP) 75, 75 (1978). 
[10] V. A. Dzuba, V. V. Flambaum , PG. Silvestrov, and O. P. Sushkov, Phys. Scr., 31, 275 

(1985). 

[11] A. Derevianko, M. S. Safronova, and W. R. Johnson, Phys. Rev. A 60, R1741 (1999). 
[12] V. A. Dzuba, V. V. Flambaum, A. Ya. Kraftmakher, and O. P. Sushkov, Phys. Lett. A 
142, 373 (1989). 

[13] M. S. Safronova, W. R. Johnson, and A. Derevianko, Phys. Rev. A 60, 4476 (1999). 
[14] V. A. Dzuba, V. V. Flambaum, and O. P. Sushkov, J. Phys. B 17, 1953 (1984). 
[15] V. A. Dzuba, V. V. Flambaum, P. G. Silvestrov, and O. P. Sushkov, J. Phys. B 20, 
1399 (1987). 

[16] Random-phase approximation; see, e.g. M Ya. Amusia and N. A. Cherepkov, Case 

Studies in Atomic Physics, 5, 47 (1975). 
[17] J. Hoffnagle, L. Ph. Rvesch, V. L. Telegdi, and A. Weis, VIII ICAP, Abstracts, p. A33, 

Goteborg, Sweden (1982). 
[18] S. L. Gilbert, R. N. Watts, and C. E. Wieman, Phys. Rev. A 29, 137 (1984). 
[19] V. A. Dzuba, V. V. Flambaum, and O. P. Sushkov, Phys. Lett. A 140, 493 (1989). 
[20] E. Arimondo, M. Inguscio, and P. Violino, Rev. Mod. Phys. 49, 31 (1977); S. L. Gilbert, 

R. N. Watts, and C. E. Wieman, Phys. Rev. A 27, 581 (1983). 



10 



TABLES 



TABLE I. Hyperfine structure matrix elements for the 6S and 7S states of 133 Cs (MHz). 





Approximation 


A 6S 


A 7S 


/ A A 
V A GS A 7S 


/R CI M7C\ 
\OD|/i| ID) 


(6S\h\7S) 
VA 6S A 7S 


- l 


V 


= l/^HF 


mm 


1424.8 


391.5 


746.9 


746.9 











(iP\h + SV\ip) 


1712.5 


469.7 


896.9 


897.1 


2.2 x 10" 


-4 






{i)\h + 8V + 8E\i/>) 


1687.8 


466.9 


887.7 


887.6 


1.1 x 10" 


-4 




= i>Br 




1952.4 


459.5 


947.2 


947.2 











(tP\h + SV\ip) 


2302.0 


541.4 


1116.3 


1116.7 


3.5 x 10" 


-4 






{ip\h + 5V + St\if>) 


2267.6 


537.7 


1104.3 


1104.5 


1.8 x 10" 


-4 




= 1>fi? 


(i/>\h + 5V\ip) 


2308.3 


542.5 


1119.0 


1119.5 


4.4 x 10" 


-4 






{ip\h + 5V + S±\ip) 


2273.8 


538.8 


1106.9 


1107.3 


3.6 x 10" 


-4 






SDpT b 


2278.5 


540.6 


1109.8 












Experiment 


2298.2 


545.9 


1120.1 









a Brueckner orbitals with £ operator rescaled to fit the energy. 

b Single, double and partly triple excitation approximation; calculations by the Notre-Dame group, 



reference [13[ 



^Reference [EQ 
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FIGURES 
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FIG. 1. Second-order diagrams for the self-energy of the valence electron (S operator). Dashed 
line is the Coulomb interaction between core and valence electrons. Loop is the polarization of the 
atomic core which corresponds to the virtual creation of the excited electron and a hole in the core 
shells. 





FIG. 2. Bruckner-type correlation diagrams for the hfs. Cross denotes the hfs interaction. The 
S operator includes second-order diagrams (Fig.l) and higher-order diagrams as described in the 
text. 
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FIG. 3. Core polarization (RPA) diagrams for the hfs in the first and second order in Coulomb 
interaction. 
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